WHAT IS CLAIMED IS : 

1 . A method for maximum a posteriori (MAP) decoding of an input 
infomiation sequence based on a first information sequence received through a channel, 
comprising: 

iteratively generating a sequence of one or more decode results starting with an 
initial decode result; and 

outputting one of adjacent decode results as a decode of the input information 
sequence if the adjacent decode results are within a compare threshold. 

2. The method of claim 1, wherein the iteratively generating comprises: 

a. generating the initial decode result as a first decode result; 

b. generating a second decode result based on the first decode result and a model 
of the channel; 

c. comparing the first and second decode results; 

d. replacing the first decode result with the second decode result; and 

e. repeating b-d if the first and second decode resxilts are not within the compare 
threshold. 

3. The method of claim 2, wherein the generating a second decode result 
comprises searching for a second information sequence that maximizes a value of an 
auxiliary function. 

4. The method of claim 3, wherein the auxiliary function is based on^e^ 
expectation maximization (EM) algorithm. 

5. The method of claim 4, wherein the model of the channel is a Hidden 

Markov Model (HMM) having an initial state probability vector % and probability density 

matrix (PDM) of P(X,Y), where XeX, YeY and elements of P(X,Y), 
pij(X,Y) = Pr(j,X,Y I i), are conditional probability density functions of an information 
element X of the second information sequence that corresponds to a received element Y 
of the first information sequence after the HMM transfers firom a state i to a state j, the 
auxiliary function being e>^pressed as: 
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'■is 3^ 



■ —a 



^ - 



8 Q(X^X/;p) = j;^(z,X;p,Y,^)log(^(z,X^Y7)), where P is a number of 

2 

T 

9 iterations, T(z, Xj , ^\) = ^\^Y[ ^u-iu (X t , ) , T is a number of information elements 



t=i 



10 in a particular information sequence, z is a HMM state sequence ij , 71 is the probability 

11 of an initial state io, X/^ is the second information sequence, X^p is a second information 

1 2 sequence estimate corresponding to a pth iteration, and Y,^ is the first information 

13 sequence. 

1 6. The method of claim 5, wherein the auxiliary function is expanded to be: 



T n n 



2 Q(X^ ) = X S Z Y,u( K ) log (Po(X., YO) + C 

t=l i=l j=l 

3 where C does not depend on X^ and 

4 y.,j(X^p) =ai(X|:;, Y;-')py(X,p, Y.) X[,,,, . Y^,. ) 

5 where a j (X| p , Yj^ ) and pj( X^^| p , Y^^j ) are the elements of forward and backward 

6 probability vectors defined as 

7 a(X|, Y;) = 7l^P(Xi,Yi),andP(X^Y,^) = ^P(Xj,Yj)l, Tiisan 

i=I j=t 

8 initial probability vector, 1 is the colunm vector of ones. 

1 7. The method of claim 6, wherein a source of an encoded sequence is a 

2 trellis code modulator (TCM), the TCM receiving a source information sequence l[ and 

3 outputting X[ as an encoded information sequence that is transmitted, the TCM defining 

4 Xt = g^Sj where Xt and I, are the elements of X[ and for each time t, respectively, 

5 is a state of the TCM at t, and gXO is a fimction relating X,, to 1^ and Sj, the method 

6 comprising: 

7 generating, for iteration p+1, a source information sequence estimate I^p^j that 

8 corresponds to a sequence of TCM state transitions that has a longest cumulative distance 

9 L(St.i) at t = 1 or L(So), wherein a distance for each of the TCM state transitions is 

1 0 defined by L(S^ - L(Sj+,) + m( Ij^j (St+,))for the TCM state transitions at each t for 





Docket No. 1999-Uo78 16 



11 t = 1 , T and the cumulative distance is the sum of m( (SJ) for all t, m( (SJ) being 

1 2 defined as 



13 m(i.(S,))=^ ^ Y,ij(irpKPc,ij(Y.|X.(S.)),foreacht= 1,2, ...,T, where 

14 X^(S() = gt(St, It (St)), n^ is a number of states in an HMM of the channel and 

1 5 Pc,u( Yt I Xt(St)) are channel conditional probability density functions of Yt when X^(S^ is 

16 transmitted by the TCM, I^p^, being set to a sequence of 1^ for all t. 

1 8. The method of claim 7, wherein for each t = 1, 2, T, the method 

2 comprises: 

A 

3 generating m( Ij (SJ) for each possible state transition of the TCM; 

4 selecting state trajectories that correspond to largest L(St) = L(St+,) + m( Ij^i (S,+,)) 

5 for each state as survivor state trajectories; and 

6 selecting I, (S^s that correspond to the selected state trajectories as It,p+i(St). 

1 9. The method of claim 8, further comprising: 

2 a. assigning L(S j )=0 for all states at t = T; 

A 

3 b. generating m( (SJ) for all state transitions between states St and all possible 

4 states Sj+j; 

5 c. selecting state transitions between the states Sj and S^+i that have a largest 

A A 

6 L(St) = L(St+, ) + m( 1,^1 (St+i)) and It^., (Sj+j) that correspond to the selected state 

7 transitions; 

8 d. updating the survivor state trajectories at states St by adding the selected state 

9 transitions to the corresponding survivor state trajectories at state St+i; 

10 e. decrementing t by 1 ; 

11 f. repeating b-e until t = 0; and 

A 

12 g. selecting allAe It (St) that correspond to a survivor state trajectory that 

1 3 corresponding to a largest L(St) at t = 0 as lj^p^.i . 

1 10. The method of claim 6, wherein the channel is modeled as ¥j(Y \ X) = 

2 PcBc(Y I X) where is a channel state transition probability matrix and Bj(Y | X) is a 
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3 diagonal matrix of state output probabilities, the method comprising for each t = 1 , 2, 

4 T: 

5 generating y,,( i;^^ ) = a,( Y/ M K^ I I^.p ); 

6 selecting an 1. (S^) that maximizes L(S , )=L(S )+ m( I,^, (8.+,)), where m( 1^ (SJ) 

7 is defined as 

8 m( i , (S J) = ^ Yu( I^p )Pj(Y, I X,(SO), n, bemg a number of states in an HMM of 

9 the channel; 

10 selecting state transitions between states S^and St+, that corresponds to a largest 

11 L(S,) = L(S,,) + m(i,,.(S,,0);and 

12 forming survivor state trajectories by connecting selected state transitions. 

1 11. The method of claim 1 0, further comprising: 

2 selecting I J (SJ that corresponds to a survivor state trajectory at t = 0 that has the 



largest L(St) as I^p^j for each pth iteration; 
comparing 1;^^ andl^p,i;and 



5 outputting I^p^, as the second decode result if I^p andl^p^j are v^thin the 

6 compare threshold. 

1 1 2. A maximum a posteriori (MAP) decoder that decodes a transmitted 

2 information sequence using a received information sequence received through a channel, 

3 comprising: 

4 a memory; and 

5 a controller coupled to the memory, the controller iteratively generating a 

6 sequence of one or more decode results starting with an initial decode result, and 

7 outputting one of adjacent decode results as a decode of the input information sequence if 

8 the adjacent decode results are within a compare threshold. 

1 13. The decoder of claim 1 2, wherein the controller: 

2 a. generates the initial decode result as a first decode result; 

3 b. generates a secorid decode result based on the first decode result and a model 

4 of the chaimel; 
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5 c. compares the first and second decode results; 

6 d. replaces the first decode resxilt with the second decode result; and 

7 e. repeats b-d until the first and second decode result are not within the compare 

8 threshold. 

1 14. The decoder of claim 13, wherein the controller searches for information 

2 sequence that maximizes a value of an auxiliary fiinction. 

1 15. The decoder of claim 14, wherein the auxiliary fimction is based on 

2 expectation maximization (EM). 

1 16. The decoder of claim 1 5, wherein the model of the channel is a Hidden 

2 Markov Model (HMM) having an initial state probability vector n and probability density 

3 matrix (PDM) of P(X,Y), where X€X, YeY and elements of P(X,Y), 

4 Pij(X,Y) = Pr(j,X,Y I i), are conditional probability density fimctions of an information 

5 element X of the second information sequence that corresponds to a received element Y 

6 of the first information sequence after the HMM transfers from a state i to a state j, the 

7 auxiliary fimction being expressed as: 

8 Q(X[ , X ) = ^ T(z, Xjp , Y7 ) log(^(z, X[ , Yi^ ) ), where p is a number of 

z 

T 

9 iterations, ^'(z, X^ , Yj^ ) = 7C P][ p;^ (Xj , Y^ ) , T is a number of information elements 

10 in a particular information sequence, z is a HMM state sequence i J , 7t is the probability 

11 of an initial state io, X[ is the second information sequence, X^p is a second information 

12 sequence estimate corresponding to a pth iteration, and Yj^ is the first information 

1 3 sequence. 

1 17. The decoder of claim 16, wherein the auxiliary function is expanded to be: 

2 Q(X^ ) = X Z Z K ) 'og (Py(X., Yd) + C 

5 where a j (X} p , Y,^ ) and Pj( ^t+i,p » ) ^® the elements of forward and backward 

6 probability vectors defined as 
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7 a(X{, Y;) = 7inP(X-,Yi),andp(X^,j,Y7,,) = nP(Xj,Yj)l^ 

8 probability vector, 1 is the column vector of ones. 

1 1 8. The decoder of claim 1 7, wherein a source of an encoded sequence is a 

2 trellis code modulator (TCM), the TCM receiving a source information sequence and 

3 outputting as an encoded information sequence that is transmitted, the TCM defining 

4 = gj(S, where X^ and 1^ are the elements of X^ and ij for each time t, respectively, Sj 

5 is a state of the TCM at t, and gt(.) is a function relating Xj, to I, and the controller 

6 generates, for iteration p+1, an input information sequence estimate I^p^j that 

7 corresponds to a sequence of TCM state transitions that has a longest cumulative distance 

8 L(St.i) at t = 1 or L(So), wherein a distance for each of the TCM state transitions is 

9 defined by L(S^+{) = L(St+i) + m( 1^^, (St+i)) for the TCM state transitions at each t for 

10 t = 1, T and the cumulative distance is the sxmi of m(Ij (SJ) for all t, m(Ij (SJ) being 

1 1 defined as 

12 m(i,(S.))=|; I; Y^jfc)logp,jj(Y,|X.(S,)),foreacht=1.2,...,T,where 

i=l j=l 

A 

13 Xt(St) = gt(St, It (St)), n^ is a number of states in an HMM of the channel and 

14 Pc.ij(Yt I Xt(St)) are channel conditional probability density fimctions of Yt when X^(S^ is 

1 5 transmitted by the TCM, p^^i being set to a sequence of 1^ for all t. 

1 19. The decoder of claim 18, wherein for each t = 1, 2, T, the controller 

' A 

2 generating m(It (St)) for each possible state transition of the TCM, selecting state 

A 

3 trajectories that correspond to largest L(S^ = L(St+i) + ni( l^^^ (S^^J) for each state as 

A 

4 survivor state trajectories, and selecting Ij+i (St+i)s that correspond to the selected state 

5 trajectories as It+i.f^i(Si^i). 

1 20. The decoder of claim 1 9, wherein the controller: 

2 a. assigns L(S j )=0 for all states at t = T; 

3 b. generates m(l^ (S^)) for all state transitions between states St and all possible 

4 states St+i; 




in 
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5 c. selects state transitions between the states S, and S,+, that have a largest 

6 L(St) = L(St+, ) + m( Ij^i (Sj+,)) and 1^^, (8^+,) that correspond to the selected state 

7 transitions; 

8 d. updates the survivor state trajectories at states St by adding the selected state 

9 transitions to the corresponding survivor state trajectories at state S^+i; 

10 e. decrements t by 1 ; 

1 1 f. repeats b-e until t = 0; and 

A 

12 g. selects all the I, (SJ that correspond to a survivor state trajectory that 

13 corresponding to a largest L(St) at t = 0 as I^p^j . 

1 21. The decoder of claim 20, wherein the channel is modeled as P^C Y | X) = 

2 PcBc(Y I X) where P^. is a channel state transition probability matrix and B^Y | X) is a 

3 diagonal matrix of state output probabilities, for each t = 1, 2, T, the controller: 



4 generates y,^ Ij, ) = a,i Y/ 1 1| „ M Y,I, 1 1^,, „ ); 

A A A 

5 selects an I t(SJ that maximizes L(S t)==L(S t+i)+^(It+i(St+i))» where m(It(St)) is 

6 defined as 



7 m(l, (S^)) = X Yt, i( lup )Pj(Y, I Xt(St)), n, being a number of states in an HMM of 

8 the channel; 

9 selects state transitions between states St and S^+j that corresponds to a largest 

10 L(St) = L(St,i ) + m( it,, (St.,)); and 

1 1 forms survivor state trajectories by connecting selected state transitions. 

A 

1 22. The decoder of claim 21 , wherein the controller selects (St) that 

2 corresponds to a survivor state trajectory at t = 0 that has the largest L(St) as I, for each 

3 pth iteration, compares I, p and I, , and outputs I, p+i as the second decode result if I, p 

4 and.;,, a. wiflu„*e'cpmpa.e threshold. 



